Magnetic-field dependence of transport in normal and Andreev billiards: 
a classical interpretation to the averaged quantum behavior 



in 

O 

o 



Nikolaos G. Fytas and Fotis K. Diakonos 
Department of Physics, University of Athens, GR-15771 Athens, Greece 

Peter Schmelcher 
Theoretische Chemie, Ira Neuenheimer Feld 229, 
Universitat Heidelberg, 69120 Heidelberg, Germany and 
Physikalisches Institut, Philosophenweg 12, Universitat Heidelberg, 69120 Heidelberg, Germany 



Matthias Scheid, Andreas Lassl, and Klaus Richter 
Institut fur Theoretische Physik, Universitat Regensburg, 93040 Regensburg, Germany 



O 

o : 

> 
(N 

m 

o 
in 
o 

-I— > 
c3 



i 

C 

O 

o 



X 



Giorgos Fagas 

Tyndall National Institute, Lee Maltings, Prospect Row, Cork, Ireland 
(Dated: February 2, 2008) 

We perform a comparative study of the quantum and classical transport probabilities of low- 
energy quasiparticles ballistically traversing normal and Andreev two-dimensional open cavities 
with a Sinai-billiard shape. We focus on the dependence of the transport on the strength of an 
applied magnetic field B. With increasing field strength the classical dynamics changes from mixed 
to regular phase space. Averaging out the quantum fluctuations, we find an excellent agreement 
between the quantum and classical transport coefficients in the complete range of field strengths. 
This allows an overall description of the non-monotonic behavior of the average magnetoconductance 
in terms of the corresponding classical trajectories, thus, establishing a basic tool useful in the design 
and analysis of experiments. 
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I. INTRODUCTION 

Ballistic transport of particles across billiards is a field 
of major importance due to its fundamental properties 
as well asphysical applications (see for example the re- 
views In such systems, a two-dimensional 
cavity is defined by a steplike single-particle potential 
where confined particles can propagate freely between 
bounces at the billiard walls. For open systems the pos- 
sibility of particles being injected and escaping through 
holes in the boundary is also allowed. As an example, 
we consider the open geometry of the extensively studied 
Sinai billiard shown in Fig. ^ Experimental realizations 
are based on exploiting the analogy between quantum 
and wave mechanics in either microwave and acoustic 
cavities or vibrating plates Q, and on structured two- 
dimensional electron gases in artificially tailored semi- 
conductor heterostructures Q,H,Q|. I n the latter case, 
the particles are also charge carriers making these nanos- 
tructures relevant to applied electronics. 

Focussing the attention on the electronic analogues, 
more recently the possibility to couple a superconductor 
to a ballistic quantum dot has been considered both the- 
oretically 0, |g and experimentally Q , so that some part 
of the billiard boundary exerts the additional property of 
Andreev reflection Q . During this process particles with 
energies much smaller than the superconducting gap A 
are coherently scattered from the superconducting inter- 
face as Fermi sea holes back to the normal conducting 



system (and vice versa). Classically, Andreev reflection 
manifests itself by retroreflection, i.e., all velocity com- 
ponents are inverted, compared to the specular reflection 
where only the boundary normal component of the veloc- 
ity is inverted. Thus, Andreev reflected particles (holes) 
retrace their trajectories as holes (particles). If, how- 
ever, a perpendicular magnetic field is applied in addi- 
tion, such retracing no longer occurs due to the inversion 
of both the charge and the effective mass of the quasipar- 
ticle resulting in opposite bending. Typical trajectories 
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FIG. 1: The open geometry of the Sinai billiard considered 
in this study. 
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are illustrated in Fig. [2J Here, we investigate the inter- 
play between trajectory bending and Andreev reflection 
and demonstrate how such effects influence the overall 
(magneto)transport properties of Andreev billiards when 
compared to their normal counterparts. 

A unique feature of this class of (quantum) mechanical 
systems is their suitability for studying the quantum-to- 
classical correspondence. In particular, much effort has 
been devoted in revealing the quantum fingerprints of the 
classical dynamics which may be paramctrically tuned 
from regular to chaotic via, e.g., changes in the billiard- 
shape. A range of theoretical tools has been used, span- 
ning the usual analysis of classical trajectories and the 
semiclassical approximation to the models of Random 
Matrix Theory and fully quantum mechanical calcula- 
tions. The main signatures of classical integrability (or 
lack of it) on the statistics of energy levels and properties 
of the transport coefficients for closed and open systems, 
respectively, have been discussed in detail in various re- 
views 0, 0, 0, 0- Discussions on modifications owing 
to the possibility of Andreev reflection appear in more 
recent studies |H 0, IS E3, E3, E3, E3, E3, EB| , mostly fo- 
cusing on the features of the quantum mechanical level 
density. 

In a similar fashion, the aim of this paper is to deter- 
mine how far a purely classical analysis may provide qual- 
itative rationalization and quantitative predictions for 
the average quantum mechanical transport properties of 
a generic billiard such as that of Fig.^ both in the pres- 
ence or absence of Andreev reflection. Indeed, by per- 
forming exact calculations for the classical and quantum 
dynamics of low-energy quasiparticles we find that the 
classical transport probabilities of electrons and holes, if 
appropriate, are in good quantitative agreement with the 
mean value (to be defined below) of the corresponding 
quantum mechanical scattering coefficients that deter- 
mine the magnetoconductance of such systems. While 
most of the previous works considered the case of zero 
or small magnetic field (such that the classical dynamics 
is not altered), we particularly analyse the regime of fi- 
nite magnetic field strengths and show that the classical 



trajectories which depend parametrically on the applied 
magnetic field suffice to describe the overall features of 
the observed non-monotonic behavior. 

The article is structured as follows. In Sec. [HI after 
a brief discussion on the details of the studied system, 
we present precise numerical results of the magnetic-field 
dependence of the transport coefficients as determined by 
the quantum mechanical scattering matrix. In Sec. IIIII 
we first discuss the model describing the corresponding- 
classical dynamics and provide an analysis for both the 
normal and the Andreev version of the Sinai billiard in 
Sections IIII Bl and IIII CI respectively. A synopsis is given 
in Sec. El 



II. QUANTUM MECHANICAL TRANSPORT 
PROPERTIES 



We consider ballistic transport of charge carriers in the 
2D Sinai-billiard shown in Fig.^under an externally ap- 
plied magnetic field. The side length of the square cavity 
is taken L = 5W where W is the width of each of the 
leads attached to the left and right of the cavity. The 
latter define source and sinks of (quasi)particles. The 
central scattering disk possesses the radius R — W, and 
it can be either a normal or a superconducting antidot. In 
the former case the antidot represents an infinitely high 
potential barrier while in the latter case it is considered 
as an extended homogeneous superconductor character- 
ized by the property of Andreev reflection Experi- 
mentally, such antidot structures have been realized in 
periodic arrangements, thus forming superlattices 0,0- 
The boundaries of the square cavity, numbered clockwise 
by the labels 1 through 4, are always normal conducting 
potential walls of infinite height. 

In the presence of a superconductor the quantum dy- 
namics of the system can be described by the Bogoliubov- 
deGennes Hamiltonian: 



H 



Ho A 

A* -m 



(1) 




FIG. 2: Typical specular (SR) and Andreev reflection (AR) 
at the circular central "antidot" of Fig. A magnetic field is 
applied as indicated. 



where the diagonal operators determine the motion of 
particles and holes, respectively, and the off-diagonal el- 
ements take care of the coupling between particle- and 
hole-like excitations. Later, in our classical calculations 
we assume perfect Andreev reflection meaning that all 
particles that hit the normal-superconducting (NS) in- 
terface are exactly retroreflected. In order to model this 
quantum mechanically we have to consider perfect cou- 
pling between the normal-conducting region and the su- 
perconductor and simulate a bulk superconductor. To 
this end, we take its size to be much larger than the su- 
perconducting coherence length £g = Hvf/^A; vp is the 
Fermi velocity. Under these conditions, it is sufficient to 
consider a stepfunction-like behavior of the pair potential 
so that A = Ao is constant inside the superconducting 
region and zero outside. We also assume that the tern- 
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perature is sufficiently smaller than the superconducting 
critical temperature so that £g does not diverge. 

For our numerical calculations we use a discretized ver- 
sion of the Bogoliubov-deGennes Hamiltonian Q resem- 
blin g th e tight-binding approximation on a square lat- 
tice [lja] . Hence, Tt becomes a matrix where only coupling 
between neighboring lattice sites is considered. The sub- 
matrix [Ho]ij has elements ti — Ep for i = j and 7y for 
nearest neighbors i and j. The Fermi energy Ep is set to 
a value that allows six open channels in the leads. The 
pairing matrix is given by [A]« = A <5y if lattice point i 
is inside the superconducting region and it is zero other- 
wise. To reproduce the correct dispersion relation in the 
continuum limit the onsite energies and the hopping 
energies 7y have to fulfill the relation e, = Y^(ij\lij> 
where denotes a summation over nearest neighbors 
j of site i, see e.g. Ref . fl6l 

In the presence of a magnetic field the hopping ener- 
gies acquire a phase according to the Peierls substitu- 
tion 0|, 7y = — exp[27ri^4 • (ri — rj-)/<I>o]- Here, A is the 
vector potential, — fj is the vector pointing from site 
j to site i and $0 = h/e is the flux quantum. In gen- 
eral, the pair potential Ao exp [i%(v4)] is also a complex 
number. However, it can be chosen real (x = 0) if the 
vector potential A is parallel to the screening currents 
near the NS interface [13. This is achieved by choosing 
the symmetric gauge A = [(B/2)z, 0, — (B/2)x] that ac- 
counts for a homogeneous magnetic field of strength B 
in y direction, perpendicular to the two-dimensional sys- 
tem. In what follows, we define as magnetic field unit the 
value Bq = tuvf / {~q e W) for which the cyclotron radius 
is equal to W . 

The transport coefficients are calculated via a recur- 
sive decimation method, as explained in Ref. 0, This 
method enables the exact computation of the full scat- 
tering matrix S n , n ' (s, H), which yields scattering proper- 
ties of quasiparticles with energy e, incident on a phase- 
coherent structure described by a Hamiltonian H(B). 
„ 2 

, )fl > (e, H) is the outgoing flux of quasiparticles along 

channel n, arising from a unit incident flux along channel 
n' . The quantum numbers n indicating open scattering 
channels are conveniently written as n = (i, a, v), where 
i indicates the leads, a takes the discrete values e and 
h for particles and holes, respectively, and v labels the 
quantum numbers associated with the quantization of 
the wavefunction in the transverse direction. As shown 
in Ref. 0, transport properties are determined by 



which is referred to as either a reflectance R (i = j) or a 
transmittance T (i ^ j) from quasi-particles of type (3 in 
lead j to quasi-particles of type a in lead i. After normal- 
ization to unity with the number of open channels Net, 




ability, while for a 



a =/= (3, P"j P (E, H) defines the Andreev scattering prob- 
(3, it indicates a normal scattering 



FIG. 3: Magnetic field dependence of the classical (solid line) 
and quantum (dots) transmission probability for the normal 
conducting Sinai billiard of Fig. The magnetic field is in 
units of the strength Bo for which the cyclotron radius is equal 
to W. 



probability. Such normalized quantities are equivalent to 
an angle-average and can be directly compared to the 
corresponding classical probabilities. 

In the remaining of the article we focus on the low- 
energy solutions of Eq. Q with quasiparticle energy 
e = 0, which is appropriate for the model of perfect An- 
dreev reflection at the NS interface. In this case, parti- 
cle and hole coefficients coincide. Hence, we adopt the 
shorthand notation R e = R ee /N ch = R hh /N ch , R h = 
R he /N ch = R eh /N c t and T e = T ee /N ch = T hh /N ch , 
T h = T he /N ch = T eh /N ch to indicate reflection and 
transmission probabilities, respectively. 

Due to interference effects, quantum scattering co- 
efficients are rapidly oscillating functions of the Fermi 
energy. Therefore, in order to remove the quantum 
fluctuations, we perform an energy average over values 
UfW/it G [6.2,6.8], which corresponds to six open chan- 
nels in the leads, N c t = Int[fci?W/7r], for each value of the 
magnetic field. The remaining parameters of our simula- 
tions are as follows. The width W of the lead is 25a m . 
For the superconducting antidot, we define the pair po- 
tential via A = %Vf/2£s by choosing £g = 8a m so that 
the diameter 2W is approximately 6 times larger than 
the superconducting coherence length. To define the tun- 
nel barrier in the case of the normal antidot, an onsite 
potential of 100 x h 2 /2maf n is added to all lattice sites 
lying inside. Note that the mesh lattice constant a m need 
not be defined explicitly if all energies are measured in 
H 2 /2ma^ n yielding 7^ = 1 for every i, j-pair. The above 
definitions are consistent with the requirements set in 
Ref. ITH1 about lengthscales, namely, $s,Af 3> a m and 
(,s/^f > 1. Here, Xf is the Fermi wavelength. 

First we consider the Sinai billiard with a normal an- 
tidot in the center acting as a potential barrier. In this 
case coefficients with a ^ (3, i.e., involving particle-to- 
hole conversion (and vice-versa) are identically zero as 
there is no Andreev reflection at the antidot boundary. 
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FIG. 4: Panel (a) shows the transport coefficients of parti- 
cles escaping the Andreev version of the Sinai billiard of Fig.Q 
without particle-to-hole conversion. Panel (b) shows the An- 
dreev reflection and transmission probabilities. In both panels 
the solid (dashed) line is the classical result of transmission 
(reflection) and the circles (squares) show the quantum trans- 
mission (reflection) coefficients normalized by N c ^. The field 
strength is given in units of Bo. Inset: semilogarithmic blow 
up for 0.75 < B/Bo < 1.25. 

Particles (holes) can be either normally transmitted or 
reflected. In Fig. |3 the smoothed transmission is com- 
pared to the classical curve (Sec. lIIIB|) revealing the same 
qualitative features. Even more remarkably, we see a very 
good quantitative agreement between both curves with 
deviations being within the amplitude of the small oscil- 
lations. Reflection is just symmetric to the transmission, 
i.e. R e = 1 — T e , as both the classical and the quantum 
calculation respect unitarity. 

Second we consider the case where the central antidot 
becomes superconducting. Andreev reflection now gives 
rise to non-zero Rh and Th coefficients as shown in panel 
(b) of Fig. ^] Upon comparison of the quantum results 
with the classical curves, we see again that they agree 
very nicely. The results are summarized in Fig. 21 with 
vertical lines indicating two distinct values of the mag- 
netic field, B c ,i and £? c .2, that are related to different 
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FIG. 5: Magnetoconductance for the normal (dots) and the 
Andreev (solid line) Sinai-shaped billiard of Fig. Q (in units 
of the number of open channels, A c h, times the conductance 
quantum Go = 2e 2 /h). The field strength is in units of Bq. 



qualitative features in the classical dynamics. The largest 
differences occur for the particle-to-hole coefficients at 
the first critical field B cl . There, the classical transmis- 
sion and reflection vanish abruptly, whereas the averaged 
quantum mechanical coefficients decay exponentially (see 
inset of Fig. HJ. However, we leave the analysis of such 
effects as well as the overall non-monotonic behavior with 
respect to the magnetic field for the next section for a dis- 
cussion under the prism of the properties of the classical 
trajectories. 

To conclude this section, we would like to show how 
the conductance, as an experimentally accessible quan- 
tity, changes when the antidot is made superconducting. 
In Fig. |SJ the magnetoconductance of a normal (dots) 
and for a superconducting (solid line) antidot is plotted. 
In the normal case the linear-response low-temperature 
conductance is simply proportional to the transmission 
T e , according to Landauer's formula Gn = (2e 2 //i) T e . 
Lambert et al |19| have worked out generalizations for 
systems including superconducting islands or leads. For 
the Andreev version of the Sinai billiard system of Fig.^ 
the conductance is given by Gs — (2e 2 /h) (T e + Rh). 
Overall, we see that in the presence of Andreev reflection 
the conductance of the system is larger than in the nor- 
mal conducting case for magnetic fields B < B c ,\. For 
larger fields the particle-to-hole coefficients vanish and 
the conductances for both cases almost coincide. 



III. CLASSICAL DYNAMICS 

A. General features 

In this section we study the classical dynamics of the 
incoming particles (we focus on electrons but the same 
arguments apply to incoming holes) for each of the two 
antidot structures described above. 
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The general form of the Hamiltonian describing the 
dynamics of charged particles inside the cavity reads 



H 



1 



2m', 



ip-q a Af 



(2) 



The index a is used to describe the possibility that the 
propagating particles are either electrons (e) or holes 
(h). This generalization is necessary for a correct de- 
scription of the dynamics in the setup with the super- 
conducting antidot. The canonical momentum vector is 
p = (j>x)Pz) = m* a v + q a A where v is the mechanical ve- 
locity, the corresponding position vector being r = (x, z). 
Charge conservation yields m* h — — m* for the effective 
masses and qh — —q e for the electric charge. The main 
property which distinguishes the two cases, i.e., nor- 
mal/superconducting antidot, is the interaction of the 
charged particle with the scattering disk. The latter 
is captured by the elementary processes illustrated in 
Fig. [21 namely, specular reflection (SR) versus the An- 
dreev reflection (AR). 

In what follows, we calculate the electronic transport 
properties by analyzing the ballistic propagation and es- 
cape of classical particles injected into the billiard via 
the opening pipe-like channels (see Fig. The initial 
conditions for incoming electrons are determined by the 
phase-space density 



Po(x,Z,V Xl V z ) 



1 r/ L ^ 

S(x + — )x 



W W 

e(z + T )-e(z-- 



8{m* e {v — vf)) cos#, 



(3) 



where G [— ■§,■§] is the angle of the initial electron 
momentum with the a;-axis and vf = ^/2Ef/iti* and the 
coordinate origin is assumed at the center of the cavity. 

The trajectories of the charged particles in the bil- 
liard consist of segments of circles with cy clotron ra- 
dius r = m* a v j \—q e B) (with v = \J v^. +v%). At non- 
vanishing external field the classical dynamics of both the 
normal and Andreev billiards is characterized by a mixed 
phase space of co-existing regular and chaotic regions. At 
B = the superconducting antidot leads to an integrable 
dynamics, (since trajectories are precisely retraced after 
retroreflection) , while the corresponding normal device 
possesses a mixed phase space. It is convenient to write 
the dynamics (collisions with the walls and the antidot) 
explicitly in the form of a discrete map. As the magni- 
tude of the velocity remains constant in time, a simple 
parameterization of the dynamics is given by determin- 
ing the position (x n ,z n ) of the n-th collision with the 
boundary and the angle 8 n of the velocity vector with 
respect to the normal of the boundary at the collision 
point taken after the collision. Here, the term boundary 
refers to the walls 1 throughout 4 and the circumference 
of the antidot (see Fig. 1). 

There are three families of periodic orbits each forming 
a continuous set that occur in the classical dynamics and 



phase space of the closed system [2(1 l2lt |2^. i.e. with- 
out leads, leaving their fingerprints in the open system 
with the attached leads. We will briefly discuss these 
periodic orbits in the following. At zero field there are 
orbits bouncing between two opposite walls with veloc- 
ities parallel to the normal of the corresponding walls. 
At finite but weak £?-field strength the periodic orbits 
form a rosette and incorporate collisions with the anti- 
dot and the walls. These periodic orbits are typical, i.e. 
dominant up to a critical field value B c 2 - For magnetic 
fields above B c 2 the cyclotron radius is so small that 
no collisions with the antidot can occur and skipping or- 
bits, describing the hopping of the electrons along the 
billiard walls, become dominant. All periodic orbits pos- 
sess an eigenvalue one of their stability matrix [23^ and 
all periodic orbits possess unstable directions. We remark 
that the above-discussed periodic orbits of the closed bil- 
liard are not trajectories emerging from and ending in 
the leads of the open billiard. However, trajectories of 
particles coupled to the leads (i.e., injected and trans- 
mitted/reflected) can come close to the periodic orbits of 
the open billiard thereby tracing their properties. This 
way the presence of the periodic orbits reflects itself in 
the transport properties. 



B. Sinai billiard with normal antidot 

First we consider the transport of electrons through 
the Sinai billiard (Fig. ^) with a normal antidot. The 
relevant quantities determining the current flow through 
the device are the transmission T e and reflection R e co- 
efficients for electrons defined as the percentage of the 
initial electrons leaving the device from the right and left 
lead, respectively. Additional quantities that are help- 
ful for an understanding and analysis of the system dy- 
namics are the mean number of collisions per incoming 
electron with the walls (1-4), (ri) w , and with the anti- 
dot, (n) a . We calculate these quantities by numerical 
simulation for different values of the external magnetic 
field B. It is convenient to use a dimensionless form of 
the classical equations of motion by employing the scal- 
ing x = £ X W and z = £ Z W for the spatial coordinates 
and t = t/lu (with uj = B$/m* e ) for the time coordinate. 
The above quantities are calculated for 100 values of the 
magnetic field strength varying from 0.02 to 2 using an 
ensemble of 10 6 different initial conditions distributed ac- 
cording to Eq. (0 for each B-field value. The magnetic 
field dependence of the coefficients T e and R e is shown in 
Fig. E^. The obtained curves are quite irregular, possi- 
bly indicating the presence of fractal fluctuations in the 
magnetoconductance of the system [24| • 

In Fig. we present the parametric dependence of 
the mean quantities (n) w and (n) a . Interestingly, the 
mean number of collisions with the walls remains con- 
stant (n) w ~ 10 for almost all values of the field strength. 
This means naturally that also the accumulated number 
of collisions of all injected trajectories with the walls is 
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FIG. 6: (a) Classical electron transmission, T e , and reflection, 
R e , coefficients for the normal antidot device of Fig. as 
a function of the applied magnetic field B. (b) The mean 
number of collisions with the boundary of the square cavity 
(walls 1-4), {n) w , and with the circumference of the antidot, 
(n) a , as a function of B. The field strength is given in units 
of B . 



independent of the field strength for the whole regime 
considered. This number is obtained by integrating the 
occupancy of the trajectories in phase space, i.e., their 
measure, over all possible velocities and the boundary of 
the cavity (defined by the walls 1-4 including the leads). 
Its invariance with respect to the field strength is of com- 
bined geometrical and dynamical origin and can be un- 
derstood as follows. The escape probability is well ap- 
proximated by P = fJ, csc / yUT,b where yU eS c is the measure 
of phase space points on the left and right leads visited 
by the escaping electrons, while [ir.b is the total measure 
involved in the dynamics of the system along the bound- 
ary defined by the walls 1-4 (including the leads). The 
corresponding integrals can be estimated as /i eS c = 2cjW 
and \iT.b — 4lCi,L. Here, q is the mean phase space den- 
sity on the leads and Cb is the mean phase space density 
on the entire boundary, both integrated over the mo- 
menta. Due to the symmetric setup of the leads relative 
to both the x and z axis, we have q = q, = c. Thus, 
P = 1/10 and the mean number of collisions with the 
wall is (n) w = l/P = 10. 

The behavior of (n) a is more complicated because 
the dynamical occupation of the antidot 's circumference 
strongly depends on the value of the external field. One 
can clearly distinguish three regimes: (i) the low field re- 
gion ranging from B w to B — B c \ w 0.55, (ii) the 
intermediate field region with B c \ < B < £? c .2(~ 1.33), 
and (iii) the high field region with B > B c ^- All three 



regions are characterized by different properties of the 
corresponding phase space. These are revealed by the 
study of the phase space structure using Poincare sur- 
faces of section (PSOS) for different values of the applied 
field B. We employ (x,v x ) sections defined by the con- 
dition z = 0. It turns out that all calculated surfaces of 
section reveal a mixed phase space. To further quantify 
our analysis we calculate the relative weight w c of those 
trajectories on the PSOS that exhibit collisions with the 
antidot. Since collisions with the antidot are the only 
possibility to obtain dynamics that is sensitive with re- 
spect to the initial conditions, w c is also a measure for 
chaoticity in phase space (25|. We partition the ener- 
getically allowed phase space on the (x,v x ) plane into 
N = 10 4 cells of equal size and define on each cell Ci the 
characteristic function hc^XjVx) as: 



1 if an orbit exists with (x,v x ) G C,; 

that hits the antidot, 
if (x,v x ) £ Ci holds for all 

trajectories hitting the antidot. 



(4) 



We subsequently approximate w c ~ (1/N)'^2 i _ 1 iid- 

The function w c (B) is shown in Fig. [7| The initial 
plateau at w c ~ 0.5 shows clearly that the system is to 
a large portion chaotic for low magnetic fields. This ex- 
plains the fact that in this range of fields the mean num- 
ber of collisions with the antidot (n) a is almost constant. 
The degree of chaos in the phase space of the system is 
large enough thereby ensuring that, with the exception 
of trajectories of negligible measure, each trajectory hits 
the circumference of the antidot. Hence, as evaluated 
in a similar fashion to (J,T,bi the total measure of phase 
space points nc,a involving the circumference is equal to 
2ircW. Following the arguments given above for (n) w , we 
estimate (n) a as (n} a /{n) w — nc,a 

yielding (n) a w 7r, which is in very good agreement with 
Fig. UJ>. In the intermediate field region the weight of 
the chaotic trajectories decreases almost linearly and be- 
comes vanishingly small in the high field region. The 
linear decrease of w c leads to a linear decrease of (n) a 
for this range of magnetic fields. Above B C 2 the cy- 
clotron radius of the electron trajectories is so small 
(r < L/4 — W/2) that no collision with the antidot is 
possible, i.e., (n) a = 0. 

We can now understand the non-monotonic behavior 
of the functions T e (B), R e (B) by considering the repre- 
sentative trajectory dynamics for various £>-values. Ac- 
cording to Fig. the low field region possesses two sub- 
regions: (i) < B < 0.4 and (ii) 0.4 < B < 0.55. 
Similarly the intermediate field region can be divided 
to: (i) 0.55 < B < 0.75 and (ii) 0.75 < B < 1.33. 
Note that the points defining the magnetic field win- 
dows, B w 0.4,0.55,0.75 and 1.33, also mark qualitative 
changes in the functions (n) a (B) and w c (B). 

At very low fields, B — > 0, we observe that R e is slightly 
larger than T e . This owes to many trajectories exhibiting 
only one collision with the antidot and reflected directly 
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FIG. 7: Relative weight of the chaotic part w c of the phase 
space as a function of the applied field B (in units of Bo) for 
the normal conducting Sinai-shaped billiard of Fig. Q 



back to the lead from where they came. The typical con- 
figuration consists of an incoming electron moving almost 
as a free particle, hitting the antidot once, and escap- 
ing from the billiard to the left lead (electron reflection) . 
Otherwise, in the low field region (i) the main process 
is the transmission of electrons. As the magnetic field 
increases, the reflection angle at the circumference of the 
antidot increases too and the incoming electron, after hit- 
ting the antidot, suffers two or more collisions with some 
of the walls 1,2,3 or 4 before escaping to the right lead 
(electron transmission). This mechanism, along with a 
significant amount of electrons that are initially emitted 
with a larger angle {\9\ > tan _1 (2/5)), hitting directly 
the upper or lower wall of the cavity (or even its right 
wall) suffering specular reflection and exiting to the right 
lead of the device, establishes electron-transmission as 
the main process in the low field region. However, it 
is evident from Fig. that the difference between T e 
and R e is quite small. There are complex trajectories 
with more than 10 collisions with the walls and 5-8 col- 
lisions with the antidot possessing a finite measure in 
phase space. These give a non- vanishing contribution 
to R e thereby maintaining its mean value around 0.45. 
In fact T e and R e fluctuate insignificantly around their 
mean values (0.55 for T e and 0.45 for R e ). As the fields 
increases above B = 0.4, trajectories with a larger num- 
ber of collisions with the walls may become statistically 
more important but as long as B < 0.55 the trajectories 
with several collisions with the antidot have still a signif- 
icant measure (see Figs. [7| and Ep, respectively). Overall, 
for 0.4 < B < 0.55, the main process is the reflection 
of electrons yielding a large difference between T e and 
R e . Typical trajectories have one or two collisions with 
the walls and a single collision with the antidot. The 



incoming electron hits the antidot once, is specularly re- 
flected and escapes from the billiard to the left lead, after 
suffering one more collision with the wall 4. 

At intermediate fields, we observe an almost monotonic 
decrease of R e and an increase of T e owing to the com- 
bined decrease of (n) a (B) and w c (B). The small plateau 
feature in Fig. \7\ around w c = 0.6 is also reflected in 
the change of slope in the transport probabilities within 
subregion (i). At its upper limit, B = 0.75, the two co- 
efficients become equal, T e — R e . Most trajectories have 
5 — 8 collisions with the walls and 1 — 4 collisions with 
the antidot. In window (ii) of the intermediate-field re- 
gion with B > 0.75, the process of transmission is much 
stronger than the process of reflection. Most trajecto- 
ries have a few (« 5) collisions with the walls and no 
collision with the antidot. A typical trajectory of the in- 
coming electron, due to the small cyclotron radius, suffers 
one collision with wall 1, three collisions with wall 4, one 
collision with wall 3 (5 collisions with the walls in total) 
and then escapes from the cavity through the right lead 
(electron-transmission). There is also the case in which 
the incoming electron misses the right lead of the device, 
and after suffering many collisions with the walls finally 
escapes to the left lead, contributing to the process of 
reflection. The same scenario is valid also for the high 
field region. 



C. Sinai billiard with superconducting antidot 

Compared to the dynamics of the cavity with the nor- 
mal antidot, the Andreev Sinai billiard with the central 
superconducting disc exhibits basic differences due to the 
occurrence of trajectories which suffer Andreev reflection, 
instead of specular reflection, at the circumference of the 
antidot. First, the complete description of the transport 
properties of the system requires the introduction of two 
additional coefficients describing electrons that escape as 
holes either to the left or the right lead; Rh (reflection) 
and Th (transmission), respectively. 

In Figs, m, and [Hp, all probabilities are plotted as a 
function of the applied field B. In Fig. |SJ;, we present 
the quantities (n) w (B) and (n) a (B) following the defi- 
nition of Sec. IIII Bl T a and R a (with a = e,h) exhibit 
irregular fluctuations as a function of B, similar to those 
obtained for the normal billiard. The function {n) w is al- 
most identical to the corresponding function obtained for 
the normal case. Qualitatively, the mean number (n) a of 
collisions with the superconducting antidot is also simi- 
lar to that in Fig. [Bp. However, while B C 2 is remaining 
the same, the value of the critical field B c \ is shifted to 
the larger value B c \ « 0.95. The former should be ex- 
pected since it does not involve any collisions with the 
superconducting disc. 

The difference in B c \ is explained by calculating the 
relative weight of the chaotic trajectories as in Sec. IIII Bl 
The result is shown in Fig. For magnetic fields up to 
B k, 0.95 a large part of the phase space of the system is 
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FIG. 8: Classical results for the billiard of Fig. Q with a 
superconducting disc in the center, as a function of the ap- 
plied magnetic field B (in units of -Bo), (a) Classical electron 
transmission, T e , and reflection (R e ) coefficients, (b) Clas- 
sical electron-to-hole (Andreev) transmission, Th, and reflec- 
tion, Rh, probabilities, (c) mean number of collisions, (n) w , 
with the boundary of the square cavity (walls 1-4) and with 
the circumference of the antidot, (n) a . 



chaotic ensuring the equal mean phase space density on 
the boundary of the square cavity and the circumference 
of the antidot. For B > 0.95 an almost linear decrease of 
w c leads to a corresponding linear decrease of (n) a . At 
B m 1.33 the chaotic part of the phase space vanishes 
due to the fact that no collisions with the defocusing 
perimeter of the antidot are possible. An additional pe- 
culiarity of the superconducting device appears for the 
intermediate field region, 0.95 < B < 1.33: All possible 
trajectories possess an even number of collisions with the 
antidot, yielding a vanishing transmission and reflection 
of holes, T/j = Rh = 0. This interesting feature is re- 
lated to the generic properties of Andreev reflection and 
is analyzed in Appendix A. 

Let us now consider the transport coefficients in more 
detail. Following the qualitative features in Figs. IHt 
and[SjD, we divide the low field region into four windows: 
(i) < B < 0.3, (ii) 0.3 < B < 0.5, (iii) 0.5 < B < 0.7 
and (iv) 0.7 < B < 0.95. Note that different regimes 
roughly coincide with different qualitative features of the 
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FIG. 9: Relative weight of the chaotic part of the phase space 
w c as a function of the applied field B (in units of Bo) for the 
Andreev version of the billiard in Fig. 



functions (ri) a (B) and w c (B), as indicated previously. 

In region (i), due to the large cyclotron radius (al- 
most vanishing curvature) most trajectories have a single 
collision with the antidot. In fact, at very small fields, 
B ~ 0.01, there is practically no collision with the walls. 
The typical process consists of an incoming electron mov- 
ing almost on a straight line hitting the antidot once, be- 
ing converted into a hole which nearly retraces the elec- 
tron path moving towards the left lead and finally leaving 
the device. Therefore, Rh in this region is larger than 
Th, T e and R e . There is however a significant amount 
of electrons emitted initially with a large enough angle, 
\6\ > tan _1 (2/5), which hit directly the upper or lower 
wall of the device suffering normal reflection and then ex- 
iting to the right lead. This process gives a finite electron 
transmission coefficient, yielding a significant contribu- 
tion to T e . For even larger emission angles, there is a 
second but less significant set of trajectories that exhibit 
specular reflections at the walls 2, 3, and 4. Particles 
along these paths escape finally from the cavity through 
the left opening, overall leading to very small values of 
R e . With increasing B, the curvature of the trajectories 
also increases so that an additional collision with the wall 
takes place. The incoming electron hits again the anti- 
dot, becoming a hole and as the curvature is increased 
the hole cannot escape from the narrow left lead, sub- 
sequently hitting wall 1, being specularly reflected and 
leaving the device from the right lead. Therefore, Th is 
increased at the cost of Rh- 

In region (ii), the curvature of the trajectories increases 
further. A typical trajectory for an incoming electron, 
after being Andreev reflected at the antidot, hits the wall 
4 and escapes to the left lead after specular reflection. As 
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a result Rh increases in this region. We also observe a 
decrease of T e owing to the increased curvature of the 
trajectories making it difficult for the incoming electrons 
to avoid the collision with the antidot and it is therefore 
harder to encounter the outgoing right lead. In region 
(iii), there is no significant variation of the observables 
T e< h and R c m- 

As the intermediate field regime defined by the lin- 
ear decrease in (n) a is approached, in (iv) we observe a 
sudden decrease of both T/j and Rh with a simultaneous 
increase of, predominantly, T e . Although the measure 
w c of trajectories with several collisions with the antidot 
is non-zero, thereby, contributing to the chaotic part of 
the phase space, the predominant part of the trajecto- 
ries do not experience Andreev reflection. In the regime 
with B c .i ~ 0.95, we encounter Th = Rh — due to 
the fact the every trajectory has an even number of col- 
lisions with the antidot (see Appendix A). The predom- 
inant part of the orbits have 5 to 7 collisions with the 
walls 1,3 and 4 thereby hopping along the boundary of 
the square cavity. Finally, the high-field region is char- 
acterized by Th = Rh = 0, and (n) a — as in the case of 
the normal antidot. 



IV. SUMMARY AND OUTLOOK 

Performing simulations of the classical and quantum 
dynamics of low-energy quasiparticles, we showed that a 
purely classical analysis may be used as an interpreta- 
tion tool for the average transport properties of generic 
normal and Andreev billiards. In particular, the para- 
metric dependence on the strength of a perpendicular 
magnetic field B was studied. As the strength increases, 
this dependence of the classical trajectories on the ap- 
plied magnetic field drives the classical dynamics from 
mixed to regular for both types of billiards. The latter 
is grossly reflected in the non-monotonic behavior of the 
magnetoconductance at intermediate fields. 

Owing to the increasing trajectory bending, a slight 
increase of the conductance Gat at small fields of the 
normal billiard is followed by a significant valley whose 
minimum defines the passage to intermediate strengths. 
Scattering of particles with the Sinai-billiard disc starts 
reducing, also triggering the relative weight w c of the 
chaotic part of phase space to shrink. Around Bq , which 
corresponds to a cyclotron radius equal the disc radius, 
skipping orbits settle in and transport properties con- 
verge towards the high-B limit. 

Turning on the superconductivity at the Sinai-billiard 
disc results in the interplay of the bending of the trajec- 
tories and the occurring particle-to-hole conversion. The 
magnetic field drives the integrable correlated motion of 
particles and holes into a mixed dynamics regime evi- 
dent by an initial tendency of Gs(B) towards Gn{B) ~ 
0.5 x (2e 2 /h) at small fields, typical for systems with 
phase space having a relative large chaotic part. Com- 
pared to the normal case, increasing B Andreev reflection 



counteracts to the reduction of w c , which occurs eventu- 
ally but at higher field strengths. Hence, realizing the 
Andreev billiard leads to qualitatively different behavior. 
Rather than a magnetoconductance dip, we note a reen- 
trance effect of the conductance towards its initial higher 
value. 

Our classical calculations provide not only a qualita- 
tive rationalization of the observed properties of the ex- 
act quantum mechanical scattering coefficients and of the 
corresponding magnetoconductance spectrum but also 
allow us to make quantitative predictions, as evidenced 
by the remarkable agreement between the classical and 
quantum values. This is ideal for the designing of exper- 
imental setups and a simple analysis of the results, since 
classical simulations are much less time consuming. 

In the present paper we studied the (energy) averaged 
transport properties, i.e., removing the quantum fluctua- 
tions. Yet the averaged quantum results contain weak lo- 
calization effects at zero and small magnetic fields. These 
quantum corrections to the averaged transmission are of 
order one (more precisely -1/4 for chaotic ballistic sys- 
tems) compared to the classical contribution which is 
proportional to N^. More specifically, for the transmis- 
sion per channel for the normal conducting Sinai billiard 
with N c h — 4 the negative quantum correction at zero 
field is expected to be below 0.1, in line with the nu- 
merical results depicted in Fig. |21 It is indeed possible 
to extend the existing semiclassical theory for ballistic 
weak localization 0, |26| to averaged quantum transport 
through Andreev billiards |27| . However, our focus was 
on the finite S-field range where weak localization ef- 
fects do not exist. Rather, conductance fluctuations in 
this regime encode additional quantum information, and 
previous results in closed and open billiards (either with 
a normal or a superconducting lead) indicate that such 
fluctuations are interwined with the underlying classical 
properties. Therefore, we envisage that our study could 
be further developed and utilized both theoretically and 
experimentally in future investigations. 
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APPENDIX A: DERIVATION OF THE FIRST 
CRITICAL FIELD 

Numerical results show that for magnetic fields B > 
B c i there are only trajectories with an even number of 
collisions with the superconducting antidot, yielding van- 
ishing transport coefficients Th and Rh- In order to derive 
the first critical field it is helpful to consider the mapping 
of the guiding centers of the trajectory arcs in the pres- 
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FIG. 10: Andreev reflected orbit. The mapping from one 
center to the following one is a point reflection. 




FIG. 11: The shaded area shows all possible locations for 
centers of orbits that connect the outer walls with the antidot. 



ence of Andreev reflection. Fig. ^3 shows two segments 
of an orbit right before and right after a collision with 
the NS-interface. The center of the second arc can be 
constructed from the center of the first one via point re- 
flection at the collision point at the NS-interface. 

Consider one segment of an orbit with radius r reaching 
from the outer wall to the antidot. The distance between 
the center of this arc and the outer wall has to be less 
than r and the same holds for the distance between the 
center and the antidot. This means that the center has 
to be inside the shaded area shown in Fig. ^] Consider 
an orbit that hits the superconductor only once. Such an 
orbit has to reach from the outer wall to the antidot and 
back to the wall. So the centers before and after the An- 
dreev reflection have to be located inside the shaded area. 
This condition is easiest to fulfill for an orbit that has its 
center at one midpoint of the inner quadratic boundary, 
as shown in Fig. 1111 The question is if for a certain radius 
the center after the Andreev reflection can no longer be 
mapped into the shaded region. 

The distance d between the antidot and the inner 
square is d = L/2 — r — Rq, where L is the side length 
of the square cavity and Rq is the radius of the anti- 



dot. Applying the cosine-theorem to the triangle shown 
in Fig.[niwe get R 2 = r 2 + (R + d) 2 - 2r(R + d) cos <j). 
So the angle <f> can be written as 



cos ( 



L 2 - AR 2 



8r 2 - ALr 



(Al) 



4r(L - 2R) 

The distance between the final point and the central hor- 




FIG. 12: M is the center of the first arc before the Andreev 
reflection. After one Andreev reflection it is mapped to M' , 
after another one it is mapped to M" , which has the same 
distance from the center of the antidot as M. 



izontal line is y = 2rsin0 = 2r-y/l — cos 2 4>- In order to 
get only one single Andreev reflection, the final point has 
to be inside the shaded region, therefore we have to claim 
y > Rq + d, which means 



2ry/T 



L 



+ r > 0. 



(A2) 



Solving this inequality for L = 5W and Ro = W we find 
a critical radius r c = 1.0505VF, which corresponds via 
r/W = Bq/B to a critical field of 



0.9519S . 



(A3) 



Up to now we have only shown that for a magnetic 
field B > B c ,i the particles hit the superconductor at 
least twice consecutively. But now it is easy to see that 
the number of Andreev reflections is indeed even. Two 
Andreev reflections in a row correspond to a rotation of 
the center around the center of the antidot, as shown in 
Fig. ^] After an even number of Andreev reflections the 
center of the arc is always located on the outer dashed 
circle, which has a radius greater than Rq + d. After an 
odd number of reflections at the superconductor it is on 
the inner dashed circle with a radius smaller than Rq + d. 
Therefore the particle can only 'escape' the superconduc- 
tor after an even number of collisions, which explains the 
fact. 
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